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Abstract. We present measurements of fundamental astrophysical properties of 
nearby, low-mass, K- and M-dwarfs from our DISCOS survey (DIameterS of COol 
Stars). The principal goal of our study is the determination of linear radii and effective 
temperatures for these stars. We calculate their radii from angular diameter measure- 
ments using the CHARA Array and Hipparcos distances. Combined with bolometric 
flux measurements based on literature photometry, we use our angular diameter results 
to calculate their effective surface temperatures. We present preliminary results estab- 
lished on an assortment of empirical relations to the stellar effective temperature and 
radius that are based upon these measurements. We elaborate on the discrepancy seen 
between theoretical and observed stellar radii, previously claimed to be related to stel- 
lar activity and/or metallicity. Our preliminary conclusion, however, is that convection 
plays a larger role in the determination of radii of these late-type stars. Understand- 
ing the source of the radius disagreement is likely to impact other areas of study for 
low-mass stars, such as the detection and characterization of extrasolar planets in the 
habitable zones. 



1. Introduction 



Advances in high-resolution astronomical techniques, particularly in long-baseline op- 
tical interferometry (LBOI), have enabled us to resolve directly the disks of nearby 
stars. Interferometric observations of nearby, low-mass, stars provide us with direct 
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stellar size and effective temperature measurements. These measurements provide 
a crucial resource in constraining stellar model atmospheres and stellar evolutionary 
models. 

For KM dwarfs, in particular early M, the onset of atmospheric convection is a 
phenom enon that is often imperfectly addressed in such models. As seen in the re- 
sults of iBerger et al.l (12006) . notable disagreements exist between interferometrically 



determined radii for M dwarfs and the ones calculated in low-mass stellar models 
in the sense that interferometrically obtained values for the stellar diameters are sys- 
tematically larger, by more than 10%, than those predicted from models. These re- 
sults confirm other obse r vations of larger-than-expecte d st ellar radii, such as tho se of 
iLopez-Morales & Ribasl(l2005h . lvon Braun et al.l(l2008h and lBovaiian et all (12008). and 



provide motivation for adjustments to be made to models in order to match the obser- 
vations . 

1.1. Observations 

In this work, we present observations made at the CHARA Array, an optical/infrared 
interferometer located at Mount Wilson Observatory, California. Our sample selection 
consists of ~K0-M4 dwarfs, complete out to ~ 6.5 parsecs and is limited only upon 
V-band magnitude and declination > 10 degrees. Our goal is to measure the angular 
diameters of these stars to better than 4%. Figure [T] demonstrates our data fit to a limb- 
darkened diameter with 1% precision to the fit. 




Figure 1. Example limb-darkened angular diameter fit to CHARA visibilities 
measured in /T'-band for GJ 892. 

Thus far, we have measured the diameters of 23 cool dwarfs, 9 K-stars and 14 
M-stars, to an average precision of ~l.5%. This more than doubles the current count of 
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low mass stars with published diameters in the literature measured with interferometry 
to better than 5% precision (see Table dB 

Table 1. Status of Stellar Angular Diameters Measured to Better than 5% Precision 





Published 


THIS WORK 


K-stars 


8 


9 


M-stars 


9 


14 


TOTAL 


17 


23 



1.2. Fundamental Properties 

The linear radii R of these stars are quickly determined by combining the HIPPAR- 
COS parallax with the interferometric angular diameter measurement. Furthermore, 
for each target we are able to calculate the bolometric flux Fbol by performing a Spec- 
tral Energy Distribution (SED) fit to all available photometry in the literature (see Fig- 
ure|2l), and again with ///P/'A/? COS' parallax this quantity is readily converted to a mea- 
sure of the absolute stellar luminosity L. Finally, the effective temperature Teff can 
be expressed in relation to observable quantities using the Stephan-Boltzmann equa- 
tion in the form Teff = 2341(Fbol/^ld-'''^^' where Fbol is expressed in units of 
10"^ erg/s/cm^ and ^ld is in units of milli-arcseconds. Generally speaking, the errors 
in radii and Teff are dominated by the errors in the measured angular diameteiEl where 

2. Discussion 

2.1. Empirical Relations 

The inclusion of our data, which more than doubles the sample size presently available 
in the literature (Table [B, supplies us with a hearty foundation for establishing empiri- 
cally calibrated relationships to the stellar radius and temperature. This is especially the 
case when we constrain the precision of the data set to only those stars possessing have 
accurately measured properties, and in the following analysis, we choose a threshhold 
of 0-6^0 < 5%. 

First, we present relations for the effective temperature of the star. In Figure |3j we 
show a preliminary fit to a 3'^'^-order polynomial with the {V - K) and {B - V) color 
indicies, both of which give a median absolute deviation from the fit of just over 100 K. 
However, since the median error of the temperature measurements is a little more than 
half of these values, improvements to the fit are warranted. An obvious improvement 



'A total of 17 diameter s that meet this criteria, which includes 8 K-stars and 9-M-stars llLane et alfcOOll : 
Segransan et al. 2003; Pi Folco et al Jl 20041 . IB erger et aIll200g : lBovaiian etani2008l : [Kervella et aL|[2008l : 
Demorv et all200a) . 

^The errors in Fbql and Timf are negligible and/or non-influential in the overall error budget. 
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Figure 2. Example SED fit. TOP: The red markers are the input photometry with 
corresponding flux errors (y-direction) and bandwidth dimensions (x-direction). The 
black crosses indicate the predicted model photometry based upon the source flux 
distribution (blue line) and transmission of the filter at that bandpass. BOTTOM: 
The fractional residuals to the above fit. Flux errors are correctly displayed, but 
bandwidth dimensions are excluded in this plot for clarity. 

would be in the quality of the infrared photometr)0. We expect, but cannot confirm, an 
improvement with the Teff versus (V-K) relation from the Teff versus (B-V) relation 
since the {V - K) color index is less sensitive to metallicity compared to the {B - V) 
color index. Additionally, in theory, a multi-parameter Temperature:Color:Metallicity 
fit would also improve the relation. However, metallicities of these stars are very diffi- 
cult to measure, and although a lot of work is being done with the M-type star metal- 
licity calibrations, the K-dwarfs are currently being ignored. 

We lastly would like to comment on the relations that can be derived with stellar 
spectral types. We code the spectral types using a simple linear scale (starting with AO, 
Al, A2, A3 ... KO, Kl, K2, K3 . . . ^ 0, 1, 2, 3 . . . 30, 31, 32, 33 . . .). Figure H shows 
the results for the relations with temperature and radius. Surprisingly, we find that there 
is a much tighter correlation with temperature and spectral type (median deviation in 
Teff ~ 75 K) compared to the color indices mentioned in the above paragraph. This 
is likely an artifact of the poor photometry input in the fits also discussed. In the fit 
to stellar radii, we can see there is a substantially large spread in values, which in 
turn produces a fit that is good to only ~ 10%. It is unlikely that this is due to the 



'The infrared colors for most these stars are saturated in 2MASS. For example, the 2MASS A^-magnitude 
for GJ 411 is = 3.501 ± 0.352) 
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evolutionary status of the stars since the universe is not old enough for these stars to 
have evolved from the main sequence. Perhaps a better explanation for the spread in 
radii is a difference in composition and/or activity in the stars themselves. 




Figure 3. Effective temperature versus (V - K) color index (left) and {B — V) 
color index (right). Filled points are new results presented here, where open points 
are published values in the literature. The Icr errors in Teff are shown, but in the 
majority of cases it tends to be smaller than the data point. A preliminary result to a 
3''''-order polynomial is shown as a solid line. 
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Figure 4. Spectral type index code (see text) versus effective temperature (left) 
and radius (right). Symbols are the same as in Figure[3] 



2.2. Theory Versus Observation 

A reccurring trend seen in the literature is that for stars between ~ 0.6-0.35Mo (spectral 
types ~ M0-M3.5), theory under-predicts the radius of a star when compa red to obser- 
vation s, for both single and binary stars alike (see Figure |5]l. For instance. Ifierger et"aLl 
(|2006) used the CHARA Array to measure the radii of six M-type stars, showing that 
current models are under-estimating their radii by ~ 15-20%, with a possible expla- 
nation that some missing source of opacity is not being accounted for and the metal- 
licity o f a star is correl a ted to the deviation of the model predictions to observations. 
In fact, iLopez-MoralesI (l2007l) evaluate this issue in full detail confirming that results 
are consistent with ot her published observations of sing le an d binary M-type stars (see 
Lopez-MoralesI IIoOtI and references therein). iLopez-MoralesI (l2007l) find that for the 



case of binary stars, the activity indicator Lx/L^oh is highly correlated with the radius 
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offset, whereas the single stars (which are several orders of magnitude less active than 
the binary star sample) have a less strong, but apparent, correlation with the metallicity. 
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Figure 5. TOP: The Mass-Radius relationship for K an d M stars. The solid black 
line is a 5 Gyr isochrone from the BCAH98 models (Bar affe et al.lll998|) . For stars 
with mass > 0.6 M©, the dashed line indicates L„,u - Hp and the solid line indi- 
cates L„,u - \.9Hp. BOTTOM: Deviation in radius versus mass for stars with radii 
measurements better than 5%. Masses for single s tars are derived from the /T-band 
mass-luminosity relation from lDelfosse et al] ( l2000i) . and assume a 10% error. 
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With the addition of our new data, which doubles the number of M-type star mea- 
surements in this range, these conclusions look a bit different in the context of single 
stars. In Figure [6l we show the fractional radius deviation versus metallicity [Fe/H] and 
magnetic activity level (expessed as fractional emission of the X-ray to the bolometric 
luminosity of a source Lx/Lbol) for single stars with mass < O.6M0. In these plots, it is 
apparent that there is still a ~ 10% offset in observed radii versus the model predictions. 
However, the introduction of these new data presented here negates the hypothesis that 
the deviation of the radius prediction has to do with the metallicity and/or activity of 
the stars. The offsets are consistent across the full ranges of metallicity and magnetic 
activity. 
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Figure 6. Radius deviation of single stars with mass < 0.6Mq versus metallicity 
(left) and Lx/Lbol (right). The filled black symbols are new results introduced here, 
whereas the blue symbols are published in the literature. 



We would like to note that this conclusion contradicts iDemorv et al. (120091 ). who 
assert that there is no evidence for a radius offset in single stars with respect to the 
models. To briefiy elaborate, their study inspects the entire mass range from 0.9 - 
0.1 M q. However, only the early-type M-stars have shown this offset (Lopez-Moralei 
l2007l and references therein). Furthermore, a closer look at their data set reveals that 
they only studied new diameters of K-stars (and repeated 2 M-stars). These choices of 
data restricted and/or reduced the sensitivity to the detection of an offset. On the other 
hand, although we confirm that the disagreement between observations and models still 
does exist for single stars, the diagnostics explored to shed light on the source of the 
problem presently reveal nothing to first order. 



3. Summary and Conclusion 

In conclusion, we would like to highlight how a large, homogeneous, and sensitive 
survey such as this one is crucial to the studies of low-mass stars. Currently, we have 
measured the diameters of 23 K- and M-type stars, with an average precision of ~ 
1.5% using the CHARA Array, a long-baseline optical/infrared interferometer. The 
magnitude and quality of observations presented for this project doubles the number 
of low-mass star diameters in the literature, and boosts the precision and sensitivity 
to the calibration of empirical relations to fundamental properties of stars such as the 
effective temperature and radius. Our study, with a combination of the highest quality 
new and extant data, confirms the disagreement between observations and models for 
single early M-stars at the level of ~ 10%. However, the diagnostics explored to shed 
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light on the source of the problem are inconclusive. Continuing effort in interferometry 
and photometry to support precise calibrations is warranted in order to push the limits 
of our data and extend our knowledge to a large number of stars. The data in hand 
enable us to further analyze shortfalls in stellar models, where they exi st, and provide 



motiv ation for adjustments to models in order to match the observations (iChabrier et al. 



2007 ). 
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